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Abstract 

Graphene, the new wondrous material, is a perfect example of a 
two-dimensional hexagonal crystal unlike any other. Here we exhibit 
some of the characteristic directional features associated with hexag- 
onal cells, emphasising the sixfold symmetry. We depict the X-ray, 
vibrational and electronic band structures to be expected in such sys- 
tems via 2 dimensional contour plots. 
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1 Introduction 



It is little wonder that graphene has captured the interest of physicists or 
materials scientists and no surprise that it is the centrepiece of the 2010 Nobel 
Physics Prize. The reason for the great excitement lies in is its remarkable 
transport properties and amazing strength plus the fact that technology has 
advanced to the stage where graphene can be produced in relatively large 
sheets, on a variety of substrates. A very bright future is predicted for 
it. Thus it possesses very low resistivity, high opacity and can furnish the 
photonics, spintronics and electronic industry with a low cost alternative to 
silicon. 

There are numerous review articles [U El [3l HI E] on this subject and a 
veritable avalanche of research papers has appeared, describing or exploring 
different properties. These concern its electronic characteristics, its vibra- 
tional properties, its edge effects, the Dirac-like properties of quasiparticles 
at the corners of Brillouin zones [51 El El E] , anomalous Hall effects, etc. Our 
aim in this paper is very modest: we wish to illustrate its beautiful direc- 
tional properties in a manner that might benefit teachers of condensed matter 
physics: the 2D hexagonal symmetry is somewhat unfamiliar and can serve 
as a nice extension of the standard square symmetry analysis [TU] . 

The paper is set out as follows. Section 2 contains the basic notation 
and its consequences for glancing X-rays. In Section 3 we discuss the os- 
cillation modes within a hexagon, clamped at it edges. These represent the 
energy levels of a particle held in the infinite hexagonal well; the results are no 
doubt familiar to the graphene experts but are probably not widely known to 
the usual quantum mechanics practitioners, who often resort to rectangular 
boundary conditions. We also consider the modes corresponding to antin- 
odes at the edges because these are subsequently needed for discussing the 
Kronig-Penney model in Section 5. Section 4 is devoted to idealised vibra- 
tional modes (longitudinal and transverse) of the atoms themselves, assuming 
nearest neighbour interactions. Here, as in Section 3, our aim is to illustrate 
the hexagonal symmetry of the motions. In the last Section we describe the 
electronic band structure, using the approximation of delta-function walls 
between adjoining cells rather than the tight-binding approximation used in 
most expositions. In order to keep the analysis simple, nowhere do we study 
edge effects (armchair, zigzag, corrugations [TT]) or multilayer graphene; for 
such elaboration we refer the reader to research reviews. We hope that this 
paper will serve as a useful extension of the usual quantum mechanical prob- 



2 



lems encountered by college students and highlight the beautiful hexagonal 
symmetry of the principal properties of graphene. 

2 Bases and X-Ray Features 

Let us set down our notation which will be used in subsequent sections. 
The basic hexagonal cell (with side length i) is drawn in Figure 1 and our 
coordinate system is centred there. The fundamental cell vectors are 

a = (3, -v^) £/2, b = (3,v^)£/2. (1) 

Any two nearby corner points such as carbon atoms at U and V can be 
taken as centres of X-ray scattering so that all lattice points are generated 
by translation: 

u^s = Uoo + + sb = (3r + 3s — 1, a/3(s — r — l))£/2 

Vrs = voo + ra + sb = {3r + 3s + 1, V3{s - r - l))i/2, (2) 

where r and s are integers. The basic reciprocal lattice, drawn in Figure 2, 
is nothing more than a scaled version of the original lattice, rotated by 30° 
with side length 47r/3-\/3^- Its fundamental vectors are 

A = (l,-V3)27r/3£, B = (1, ^3) 27r/3£. (3) 

Then, according to the Von Laue conditions for constructive interference, 
whenever the change of wave number is {R and 5* are integers) 

Ak = RA + SB = {R+ S, {S - R)V3) 27r/3£, (4) 

X-rays directed along the graphene plane will of course be diffracted by the 
lattice. However because there are two atoms decorating each unit cell there 
is a modulating factor 1 + exp[(i? -|- 5") 27r/3] cx 2 cos[7r(i? + S)/3] apart from 
an irrelevant phase; this means that diffraction will be greatest when R + S 
is exactly divisible by 3 but that there are no missing orders of diffraction. 

3 The hexagonal well 

The literature on the equilateral triangular drum is old indeed; it goes back to 
Lame [I2] and Pockels [13] . A very nice series of reviews on this topic has been 
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written by McCartin [13], discussing the completeness, orthogonality and 
multiplicity of the solutions, subject to Dirichlet and Neumann conditions at 
the edge of the triangle. It is not hard to adapt the solutions to the case of 
the hexagon and we shall do so shortly. (As a matter of fact the solutions of 
the equilateral triangle are embedded in a one-sixth section of the hexagon.) 
Here is an outline of the main steps. 




Figure 1: Basic cell vectors. Any two adjoining C atoms such as U and V 
can be taken as the basic cell components. 



First of all introduce triangular coordinates, u, v, w centred at the origin 
with M + f + w = 0, as illustrated in Figure 3. The maximum values of these 
coordinates within the hexagon is the radius of the incircle r = £a/3/2 and 
we see that in terms of the normal Cartesian coordinates, 



u 



-y, 



V — w 



(5) 



The second step is to note that the Schrodinger or Helmholtz equation 
for a free particle of energy E in the cell, namely 



+ 



ilj{x,y)=0; 



2ME 



may be transcribed into triangular coordinates: 



n 



2 ' 



— + 3^^ + P 
du^ d{v — wY 



^ix,y) = 0, 



(6) 



cleverly allowing for cyclic sums of separable solutions of the form. 



(v — w) (v — w) 

ip{u^v,w) = [cos(A;i'u)+C sin(A;i'u)] [cos(A;2— A^s) — h-D sin(A;2— /^s) — - — ] 
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(V — Wj (V — Wj 

+ [cos{k2u) + C sm{k2u)] [005(^3- fci) ' +£> sin(A;3- fci) ' ' ] 

o o 

(v — w) (v — w) 

+ [cos(/i;3M) + Csin(/c3M)] [cos{ki-k2) ^ ' +Dsm{ki-k2) ^ ' ] (7) 

o o 

where ki + k2 + ka = and k^ = 2{kl + kl + A;|)/3 = A{kl + kl + kik2)/3. 
The third step is to pick out symmetric or antisymmetric solutions under 
reflection about a vertical diagonal; this selects out sin or cos functions in 
fact. Typically one finds [H] the (Dirichlet) quantized solutions fcj = vrrij/r = 
/VSi- the case of interest for this section: 

. vr^iM . 7i{n2 - n3)iv - w) 

WD sym/anti = Sm . COS / Sm h 

r Sr 

'7in2U nin^ — ni)(v — w) 

sm . cos / sm h 

r 3r 

. T^n^u . TT{ni - n2){v - w) 

sm . cos / sm , (81 

r 3r 

where are non-zero integers subject to rii + n2 + = 0. It suffices to take 
ni > n2 > to obtain distinct modes (in the antisymmetric configuration ni 
and n2 need to be different for a nonvanishing answer). Hence the quantized 
energy levels are given by 

^ni,n2,n3 = + ^2 + nl)/9f = IOtt^K + ^2 + n,n2)/9f. (9) 

While we are discussing this topic let us note [H] the Neumann solutions 
(effectively antinodes along the hexagon edges) as we will require them later: 

-KUiU . 7r(n2 - n^){v - w) 

IpN sym/anti = COS . COS / Sm h 

r Sr 

'7in2U 7i(n3 ~ ni)(v — w) 

cos . cos / sm h 

r 3r 

Tinsu . TT{ni - n2){v - w) . . 

cos . cos / sm . (10) 

r 3r 

We have illustrated some of the lowest few modes in Figures 4 to 6 in the 
form of contour plots rather than 3D plots as they are easier to comprehend 
(the lightest shaded ares are the maxima while the darkest correspond to 
minima). They show the lovely hexagonal features expected of the vibrations. 
The fundamental Neumann mode is depicted in Figure 7 even though it is not 
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Figure 4: Fundamental mode (1,1) of the unit hexagonal well. 

required for the present purposes. [It should be noted that triangular cuts 
within the hexagon describe the modes for an equilateral triangle.] McCartin 
has analysed the orthogonahty and completeness properties of solutions (8) 
and (10) so we refer the reader to his review article. 
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Figure 5: Antisymmetrical overtone (2,1) of the unit hexagonal well. 




-i.o -Q.S 0.0 as ID 

Figure 7: Basic Neumann mode (1,0) of the unit hexagonal well. 

4 Vibrational Features 

Let us begin by tackling vertical (out of plane) transverse oscillations. Re- 
garding the graphene sheet as some sort of trampoline under fixed tension 
and considering only nearest neighbour interactions, the downward acceler- 
ation on a C atom located at u^^ is determined by the adjoining vertical 
displacements relative to the atoms at v^^, Vj,s_i, v^.i^. Hence vibrations are 
described by the simple harmonic equation, 

(fUrs/dt^ = -n\3Urs - Vrs - " K-ls), (H) 

where f2 is a characteristic frequency associated with the "tension" of the 
links and how strongly or loosely the C atoms adhere to the substrate. Sim- 
ilarly for the twin atom located at v^^ one obtains 

Looking for plane wave solutions [k = (/e^,, ky) below], at the positions sum- 
marized by (2), 

Urs = t/exp[z(k.Urs) — ojt], Vrs = ^ exp[z(k.Vrs) — ut] 

one arrives at two coupled equations for the amplitudes of the cell atoms at 
U and V: 
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uj^V = n''[?,V -Ue-'^-^ -2Ue'^-^'\os{kylVl/2)l (13) 
producing the eigenvalue equation 

(wVSf^^ - 1)^ = (1 + 4cJ + 4c^Cj,)/9; = cos(3A;^£/2), Cy = cos{V3kyi/2) 

or dispersion relation between frequency and wavenumber for transverse ver- 
tical acoustic (-) and optical (+) modes: 

u± = n[3±^l+4c,Cy + 4cl]. (14) 

This dispersion relation for the two modes is shown in Figure 8. In the long 
wavelength limit, 

U^_/n^ C^3f{kl + kl)/A, so VacousUc = 2n/V3i, (15) 

whereas uj+ ~ ^/6Vl for optical vibrations. The acoustic and optical modes 
remain separated when the wave vector k is directed along x even at the edge 
of the Brillouin zone, but they merge when k is directed along y and a corner 
(a so-called Dirac point) of the zone is approached. It is the instabilities in 
the vertical motion that apparently tend to make graphene sheets ripple or 
even curl up into tubes. 

Next we shall study the oscillations of the C atoms within the graphene 
plane. This time we are dealing with two-dimensional vector displacements 
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U = {Ux, Uy), V = {Vx, Vy) for the atoms located at positions determined by 
(2), rather than the vertical displacements considered previously. For small 
vibrations we assume that the restoring force is proportional to the extension 
in the direction of the link; so if one defines three unit vectors ei 2,3 along 
each of the links, the total restoring force of the atom at U summed over the 
three links is 

[(Vrs - Urs).ei]ei + [(V^-is - Urs)-e2]e2 + [(V^s-i - U^s)-e3]e3 

apart from an overall factor. Resolving the ei into Cartesian components, 
we get the equation pair for the U atom {zu is the horizontal characteristic 
frequency) , 

Uxrs/'^'^^-^Uxrs/'^ + Vxrs + iVxr-ls + Vxr»-l)/4: + V3{Vyrs-l-yyr-ls)/'K'^Q) 
Uyrs/rU^^-SUyrs/2 + 3{Vyr-ls + Vyrs-l)/^ + VSiVxrs-l-Vxr-ls)/'^- (17) 

A similar equation applies to the vector displacement V of the V atom along- 
side: 

Vxrs/ 'CO^— — 3ViErs/2-|- 6^2:^8 + (f^xr+ls + f4 rs+l)/4 + "\/3(f/yrs+l — ^yr+ls)/48) 
Vyrsl^'^= -3K/rs/2 + 3(K/r+ls + ^s+l) /4 + V^(f/s r s+1 " r+1 s) /4. (19) 

For vibrations of circular frequency uj and wave number k, one arrives at the 
secular equation. 



w 





1 + Cye-^^/2 


-iV^Sye-'^/2 \ 


/ Ux 







w 


— i\/3sye~'^ 


2,Cye~'^/2 


Uy 


1 + cye^ji 


iy/^Sye'^/2 


W 









iV^Sy^^/2 







W y 


\ VyQ^^'^^ 


1 



= 0, 



(20) 

where W = / Aw^ -3 / 2 , x = 3kJ/2, Cy = cos(\/3A';,//2), Sy = sin(\/3A-,//2) 
and w stands for a characteristic frequency of planar oscillations, which is 
not necessarily equal to the value Q for vertical oscillations because of the 
the substrate's influence. 

This leads to the set of roots (dispersion relations) : 

oj/w^O, 2V3, (3-^l + 4c^c^ + 4c2), (3 + ^l + 4c^c^ + 4c2). (21) 

The first two roots may be ignored (the first is static displacement and the 
second is dispersionless). However, the third and fourth planar roots are 
strikingly similar to the vertical motion roots (14) (acoustic and optical), 
except that the oscillation frequency is not necessarily the same. 
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5 Kronig-Penney electron band model 



We now turn to a model of the electronic energy bands which is a gen- 
eralization of the one-dimensional Kronig-Penney scheme - instead of the 
tight-binding approximation which is more popular. To that end, we en- 
visage a series of 5-function potential walls lying along the carbon links as 
that is where the electron clouds are at their most prominent. Any stray 
negative charge carrier would encounter these repulsive barriers, except of 
course at the cell corners where they would feel the attraction from the C 
nuclei. Therefore as our model we may envisage loose electrons as moving 
freely within the cells until they hit a cell edge, where the wave function 
discontinuity across the cell boundary is proportional to the wave function 
there (which is itself continuous). This should mimic the periodic properties 
and symmetry features of more realistic models. 

Consider then the following parametrizations of the symmetric wave func- 
tion ^0 for the cell round the origin, drawn in Figure 1, and the wave function 
^1 for the cell below it (centred at x = 0, y = —\/3f,). One must allow for 
both Dirichlet and Neumann boundary conditions - the analogues of sm{kx) 
and cos(A;a;) that are included in the one-dimensional analysis. Here we do 
not consider antisymmetric solutions, since the model merely serves as an 
example to illustrate the hexagonal features of the bands - the true nature 
of the bands can in fact only be revealed through experiments. 

*j = [Ajsm{kiu) + BjCos{kiu)].cos{ki + 2k2){v - w)/3 + 
[Aj sin(/c2M) + Bj cos(/c2w)]. cos(2/ci + k2){v — w)/3 — 
[Aj sm{{ki+k2)u) — Bj cos {{ki+k2)u)]. cos (/ci — k2){v — w)/3,{22) 

for cells J = or 1. There are the (dis)continuity conditions which we can 
apply at V — w or X — 0, 

{d-^,/du)\^=r-{d-^o/du)\^=r = {V/r)-^o\u=r: (23) 

as well as the Bloch periodicity conditions, 

{^^^/^u)U=r = e''"i'.{d^o/du)\u=-r , (24) 

where q is the (real) quasi-momentum and V corresponds to the size of the 
5 function barrier. Wc also remind the reader that the electronic energy is 
given by E ^ {2h'^ /3m'^)[kl + /ci + kik2]. 
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The equations (23) and (24) determine the allowed energy bands from 
the reality of q via the conditions' combination: 

(1/2 + cos(2gr)) [{a + 13) cos(a - (3) - (3 cos(2a + (3) - a cos(a + 2/3)] + 

a(cos a — cos 2a) + /3(cos (3 — cos 2/3) — (a + /3) [cos(a + (3) — cos(2a + 2/3)] 

+ V [2 sin a - sin 2a + 2 sin /3 - sin 2/3 - 2 sin(a + /3) + sin(2a + 2/3)] = 0, 

where a = kir, (3 = k2r. With ki and ^2 inclined at 120° with respect to one 
another and — 1 < cos(gr) < 1, one may plot the allowed regions of fcj space 
for the energy bands. These are drawn in Figures 9 and 10, where we see that 
if we follow the k^ axis, the bands get narrower as V increases [I5j. Indeed 
as y — oo the ratio B/A ^ in equation (22) and the k values become 
quantized as in the infinite hexagonal well. More generally the energy bands 
lie between concentric circles at low energies because of the nature of the 
expression fc^ + fcg + ^3 with triangular axes. For the tight binding model we 
refer the reader to reference [lO] . 




Figure 9: Allowed band (shaded region) for V = 1. 
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Figure 10: Allowed band (shaded region) for V = 10. 

6 Conclusions 

Our aim in this paper has not been to provide an exhaustive analysis of all 
properties of graphene, not even in summarized form. Plentiful reviews exist 
which do just that. Rather our purpose has been to highlight the symmetry 
features to expect from this remarkable material both in the electronic and 
vibrational characteristics, in the idealized situation that one has an infinite 
lattice of it. In that way, the directional properties of graphene come much 
more into evidence. In practice of course we expect the symmetry to be 
ruined by boundary effects, though the departure is probably significant only 
near the zigzag and armchair edges. 

Graphene is a perfect vehicle for teaching two-dimensional aspects of solid 
state physics and makes for a very useful variation on the standard theme of 
square lattices. University teachers can use the material to illustrate many 
aspects of condensed matter and enlarge the students' horizons. 
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